Electrical transport in two dimensional electron and hole gas on Si(001)-(2x 1) surface 
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Si(001)-(2x 1) surface is one of the many two-dimensional systems of scientific and applied interest. 
It has two surface state bands (1) anti-bonding tt* band, which has acceptor states and (2) bonding 
7r band, which has donor states. Due to its asymmetric dimer reconstruction, transport through this 
surface can be considered in two distinct directions, i.e. along and perpendicular to the paired dimer 
rows. We calculate the zero bias conductance of these surface states under flat-band condition and 
find that conduction along the dimer row direction is significant due to strong orbital hybridization. 
We also find that the surface conductance is orders of magnitude higher than the bulk conductance 
close to the band edges for the unpassivated surface at room temperature. Therefore, we propose 
that the transport through these surface states may be the dominant conduction mechanisms in 
the recently reported scanning tunneling microscopy of silicon nanomembranes. We also calculate 
the zero bias conductance under flat-band condition for the weakly interacting dangling bond wires 
along and perpendicular to the dimer row direction and find similar trends. Extended Hiickel theory 
is used for the electronic structure calculations, which is benchmarked with the GW approximation 
for Si and has been successfully applied to Si systems in past. 

PACS numbers: 73.20.-r, 73.40.-c, 73.63.-b 



I. INTRODUCTION 

Two dimensional (2D) electronic systems [l[ have been 
of great interest. In metal-oxide-semiconductor (MOS) 
technology, the 2D electron and hole gas is the key for 
functioning of the modern integrated circuits. Also, in 
high-T c materials, the 2D cuprate planes are thought to 
be responsible for superconductivity. Graphene is an- 
other noteworthy example of 2D transport and has at- 
tracted tremendous interest recently. 

The electronic structure of the Si(OOl) surface has been 
a topic of study for decades. Although, some aspects of 
this particular surface are still controversial, overall it is 



TABLE I: Details of the model systems. The models consist 
of sixteen atomic layers shown in Fig. 1 to eliminate quan- 
tization effects and hence to obtain bulk Si band gap. The 
total number of Si atoms (excluding H atoms on top and back 
surface) per unit cell are also given. As shown in Fig. 1, a 
(2x1) unit cell has two repeating unit cells in the [110] direc- 
tion (perpendicular to dimer row) and one unit cells in the 
[TlO] direction (along dimer row). 



Model Unit cell surface Unit cell" Atoms 6 



I 


H:Si(001)-(2xl) 


2x1 


16 


II 


Si(001)-(2xl)-AD c 


2x1 


16 


III 


Paired dangling bond wire c 


4x1 


32 


IV 


Si(001)-(2xl)-AD d 


2x1 


16 


V 


Paired dangling bond wire d 


4x1 


32 



"Multiples of 3.84 A - lattice constant of bulk unit cell. 
'Number of Si atoms only per unit cell. 
c in [110] direction, 
''in [110] direction. 



a well understood surface [2j, |3| . An unpassivated Si sur- 
face with paired asymmetric dimer (AD) reconstruction 
introduces two surface bands which have a 2D character 
0, 0| . The band corresponding to the anti-bonding state 
(tt* band) has acceptor states [||. Similarly, the band 
corresponding to the bonding state (tt band) has donor 
states [4|. The surface state density is about 10 15 cm~ 2 
for 7T* and tt bands. Since the tt* band is fully unoc- 
cupied and the 7r band is fully occupied, they do not 
contribute to the transport in the absence of doping and 
surface band bending. Apart from this, the tt* and tt 
states are localized on the bottom and top dimer atoms, 
respectively Q. Moreover -the wavefunction is localized 
within 10 A of surface @, which is a very important 
length scale for the system under study. 

Furthermore, these two bands may start conducting 
with doping and/or surface band bending. In such a sce- 
nario, the transport through the tt* and tt bands can be 
thought of as electron and hole transport, respectively. 
Recently, Zhang et al. [6| reported scanning tunneling 
microscopy (STM) of silicon nanomembranes with about 
10 15 cto~ 3 p-type doping. With the atomic resolution 
for this system in STM, they propose that surface dop- 
ing enabled by thermal excitation of electrons in tt* band 
results in holes in the valence band, which enables the 
hole conduction inside the valence band. Zhang et al. Q 
place the tt band edge few tenths of eV below the valence 
band edge. We find that over the transverse Brillouin 
zone, the bottom of tt* band is about 0.4 eV below the 
conduction band edge (E c ) and the top of the valence 
band is about 0.2 eV above the valence band edge (E v ). 
Additionally, our model predicts that the zero bias con- 
ductance through these surface states is at least three 
orders of magnitude higher than that of the bulk for the 
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FIG. 1: (color online) Ball and stick models. Unpassivated paired asymmetric dimer (AD) surface and hydrogenated symmetric 
dimer surface. Top four layers, represented by J, are relaxed due to surface reconstruction. Bottom twelve, represented by f, 
are bulk layers. The back surface is hydrogenated to eliminate any inadvertent DB induced states. The top Si atom of the 
AD surface is shaded as black and the bottom Si atom is shaded as gray. The atoms below the first mono-layer of Si surface 
are represented by white circles. The direction parallel to the dimer row is conveniently referred to as [TlO], whereas the one 
perpendicular to the dimer row is labeled as [110]. Atomic coordinates for these structures are discussed and reported in Ref. 
Q. Atomic visualization is done using Gauss View Q. 



doping being used in STM Therefore, we propose 
that the surface state transport alone may be the domi- 
nant conductance mechanism in silicon nanomembranes. 

We calculate transmission and zero bias conductance 
under flat-band condition for these surface states in di- 
rections parallel to and perpendicular to the dimer row 
referred to as [TlO] and [110], respectively. We also calcu- 
late transmission and zero bias conductance under flat- 
band condition for paired dangling bond (DB) wires sep- 
arated by hydrogenated DB wires in above-mentioned 
directions. To the best of our knowledge, these are the 
first atomistic transport calculations for the systems un- 
der consideration. 



II. MODEL SYSTEMS 

In a previous study we report that an isolated 
unpaired DB on an otherwise perfectively hydrogenated 
Si(001) surface will only affect its neighboring Si atoms 
within 10 A. Same is the case for the paired DBs on 
the Si surface [||]. The it* and ir states, respectively are 
localized on the bottom and top Si atoms in the paired 
AD. In the [TlO] direct ion, the Si dimer atoms, on which 
7r* and 7r states are present, are only 3.84 A apart, while 
in [110] direction, about 7.68 A apart. The hybridization 
anisotropy in [TlO] and [110] directions results in different 
band- widths [3| ■ The extent of the wavefunction overlap 
also affects the transport in these directions. 
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FIG. 2: (color online) Transmission calculations for Model I, 
II and III. Hydrogenated Si(001) surface (model I) gives bulk 
Si band gap of 1.18 eV. Unpassivated Si(001) surface (model 
II) gives large transmission in the dangling bond row direction 
([TlO]) for both the 7r* and the 7r states because the dimer 
atoms on which these states are localized are 3.84 A apart. 
Transmission through the paired dangling bond (model III) 
wire along the dangling bond row direction, is about half of 
that in model II due to half of the total dangling bond states 
per unit transverse length. 

In this paper, we consider five model systems for the 
transport calculations described in Table I and Fig. 1. 
The transport direction is shown by arrow and the unit 
cell used for each model is shown by the dashed line. 
Model I is a hydrogenated surface and is expected to 
have bulk Si band gap. We calculate transport for this 
model in [TlO] and use it as a reference for the other mod- 
els. Model II and IV are unpassivated surfaces and the 
transport is calculated in the [TlO] and the [110] direc- 
tions, respectively. Model III and V are paired DB wires 
separated by hydrogenated wires and transport is calcu- 
lated in [110] and [110] directions, respectively. The size 
of the unit cell used and number of atoms per unit cell 
are also given in Table I. Since DBs interact within 10 A, 
enough neighboring unit cells are included to calculate 

Hamiltonian (H) and overlap (S) matrices for the E( k ) 
calculations. 



III. THEORETICAL CALCULATIONS 

We use extended Hiickel theory (EHT)for the elec- 
tronic structure calculations as in Ref. [3]. EHT pre- 
scribes a semi-empirical tight-binding procedure using 



FIG. 3: (color online) Transmission calculations for models IV 
and V. For model IV, transmission through Si surface states 
in the direction perpendicular to dimer row direction ([110]) 
is smaller than the one along dimer row direction (model II) 
due to the reduced hybridization between surface states - DB 
atoms being 7.68 A apart. Conclusion for the DB wire, for 
which the transmission is calculated perpendicular to dimer 
row direction (model V) is the same. 

a Slater type orbital non-orthogonal basis set. For Si, 
EHT is benchmarked 0] with GW approximation and 
thus gives correct band structure features, such as band 
offsets and dispersions. Transferable EHT parameters 
used in this paper are taken from Ref. 0] and are sum- 
marized in Table II. We use one orbital (Is) basis set for 
the H atom, whereas a nine-orbital (3s, 3p and 3d) basis 
set is used for the Si atom. 

In order to calculate the transport properties for a 
2D channel, we first calculate the E(k) diagrams over 
the transverse Brillouin zone for each wavevector in the 
transverse direction (kt) by transforming the real space 
Hamiltonian (H) and overlap (S) matrices to reciprocal 
( k ) space: 

N 

H(k) = H mn e a -^- r ^ (1) 

m=l 
N 

S(k ) = Y S mn e a -^- r ^ (2) 

m— 1 

where k = (k t ,k) - k is the wavevector in the trans- 
port direction. For each transverse wavevector (kt), the 
system thus becomes one-dimensional (ID) and hence its 
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FIG. 4: (color online) Zero bias conductance under flat-band condition. For model I (hydrogenated surface with transport in 
the dirtier row direction), the conductance decreases exponentially inside the band gap due to Fermi tail. Close to the band 
edge, the zero bias conductance is orders of magnitude higher for model II (unpassivated surface with transport in dimer row 
direction) than that of model I. The zero bias conductance conductance for model III (DB wire with transport in dimer row 
direction) follows the same trend as that of model II with some quantitative differences. The zero bias conductance for the tt* 
states of model IV (unpassivated surface with transport perpendicular to the dimer row direction) is two orders of magnitude 
higher and less than one order of magnitude higher than that of model V (DB wire with transport perpendicular to the dimer 
row direction) inside the band gap. 



transmission is independent of the dispersion. This en- 
ables transmission to be calculated numerically by count- 
ing band-crossings at a particular energy for each kt- Fi- 
nally transmission per unit length is calculated by sum- 
mation over transverse k t as follows: 



T(E) 



k t J 



dk t T{E u k t 



(3) 



Additionally, the zero bias conductance (G ) under flat- 
band condition at a finite temperature is calculated as 
(see appendix): 



G 



2q 2 1 



I, l, T J dE T ( E ) [1 + e (£-Mo)/fcT]2 



(4) 



IV. DISCUSSION OF RESULTS 



We report the transmission per unit transverse length 
through model I in Fig. 2, which shows clean bulk band 
gap of 1.18eV as expected due to the hydrogen passiva- 
tion and is used as a reference for the other calculations. 
For the AD surface (model II), transmission through the 
tt* and 7T states is shown in Fig. 2 along the [110] di- 
rection. The atoms are 3.84 A apart on which 7r* and 
7T states are localized. Since wavefunctions are extended 
up to 10 A, hybridization is strong and hence the trans- 
mission is large. The transmission through the tt* state 
starts increasing about 0.5 eV below the conduction band 
edge. The maximum transmission through the tt* state is 
about 3 x 10 9 /m, which gives a zero bias conductance of 
about 0.2 S/fim according to Eq. 4 and is shown in Fig. 
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TABLE II: EHT parameters used for H and Si atoms. 
Keht = 2.3 

Orbital E on - aU e{eV) C l C 2 6 i^ 1 ) 6 {A~ 

H:ls -16.15180 0.53558 0.88560 

Si:3s -18.10264 0.70366 1.83611 

Si:3p -11.25298 0.02770 0.98313 0.78901 1.70988 

Si:3d -5.347060 0.68383 0.46950 0.68292 1.72102 



more, we leave the role of defects, the electrostatic ef- 
fects of the dopant atoms and dephasing due to electron- 
phonon scattering for future work. Interestingly, it has 
been suggested that atomic steps on Si(lll) surface can 
decrease conduction by two orders of magnitude [l(| • 

V. CONCLUSIONS 



4. Additionally, the transmission through the 7r state 
is about 0.2 eV above the valence band edge. For the 
model III, the transmission is again large because atoms 
corresponding to tt* and it states are 3.84 A apart. But 
since the overall density of dangling bond is half as that 
of model II, the transmission is almost half for both the 
states. However, there are some differences in charac- 
teristic features due to varying band offsets at various 
transverse wavevectors (fct) for models II and III. 

Models IV and V are the same as models II and III, 
but the transmission is calculated in [110] direction. As 
shown in Fig. 3, the overall transmission is smaller even 
inside the conduction and the valence bands, because Si 
atoms in the top four surface layers are further apart due 
to asymmetric reconstruction. Apart from this, for model 

IV, since atoms corresponding to 7r* and tt states are now 
7.68 A apart, there is small hybridization between atoms 
on which 7r* and n bands are localized and hence surface 
state hopping is quenched. The maximum transmission 
inside band gap for model IV is about 0.75 x 10 9 /m, 
which gives a conductance of about 0.05 S/m. For model 

V, the unpassivated bottom and top dimer atoms in the 
transport direction are about 15.36 A apart, which re- 
sults in very small hybridization of the 7r* and the 7r 
states. Therefore, the transmission through these states 
for model V is even smaller. 

We show the calculated zero bias conductance under 
flat-band condition at room temperature on a linear and 
a logarithmic scale in Fig. 4. The conductance for model 
II and III is orders of magnitude higher than the model 
I inside the band gap. For model IV, the conductance is 
about two orders of magnitude higher for the 7r* states 
and less than one order of magnitude higher for the tt 
states than those of model V. Furthermore, the zero bias 
conductance decreases exponentially inside the band gap 
due to the Fermi tail. For the experimental conditions 
used in Ref. the band bending is small and hence 
based on our flat-band zero bias conductance calcula- 
tions, we propose that the surface state hoping is the 
dominant conduction mechanism in the recentlyreported 
STM of lightly doped silicon nanomembranes [6J] . 

Since the transport is coherent and the channel is in 
the ballistic regime, the transmission is independent of 
the channel length. Hence, the results reported provide 
an upper limit of transmission and conductance for these 
surface states. In this study, effects ofgate voltage and 
realistic contacts are ignored. In Ref. [9(, we provide the 
basic formalism to perform these calculations. Further- 



We have studied the transport through Si(001)-(2xl) 
surface states with asymmetric dimer reconstruction in 
the dimer row and its perpendicular direction. We find 
significant transmission in the dimer row direction due 
to strong hybridization between dangling bonds 3.84 A 
apart. However, in the direction perpendicular to the 
dimer row, bottom and top dimer atoms are about 7.68 
A apart, on which 7r* and 7r states are localized, their hy- 
bridization is weak. We find similar trends for the paired 
dangling bond wires separated by hydrogenated wires. 
Furthermore, the zero bias conductance under flat-band 
condition at room temperature around the band edges is 
at least three orders of magnitude higher for these sur- 
face band in comparison with that of passivated surface. 
Therefore, we propose that conduction through the sur- 
face states may be the dominant conduction mechanism 
in the recently reported scanning tunneling microscopy 
results of silicon nanomembranes on insulator. 



Appendix 

We derive Eq. 4 as follows. For a finite voltage V, the 
difference of the Fermi function is given as: 



J1 31 l + e (E-» )/kT) l + e (E-» +qV)/kT) 

e (E-n )/kT^ e qV/kT _ 1 j 
h- !l= [1 + e (E-v. )/kT^l + e (B- Mo +eV)/fcT] 

For \qV\ -C fcT, using Taylor series expansion and re- 
taining the first term for e qV ^ hT , it can be shown that: 

qV e {E-^)/kT 
h ^ ~ kT[l + e (£-M/fcT]2 

Starting with Landaiier's approach, 

I=^f dET{E)(h- f 2 ) 

it can be shown that zero bias conductance at a finite 
temperature is given as: 

_ I _ 2q 2 1 f e (E-^)/kT 

G °-y- ~kf J dE T[E) [l + e(s-Mo)AT ]2 
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